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Abstract. The Jordan-Wigner transformation is applied to study magnetic properties of the quantum 
spin-| XX model on the diamond chain. Generally, the Hamiltonian of this quantum spin system can be 
represented in terms of spinless fermions in the presence of a gauge field and different gauge-invariant ways 
of assigning the spin-fermion transformation are considered. Additionally, we analyze general properties of 
a free-fermion chain, where all gauge terms are neglected and discuss their relevance for the quantum spin 
system. A consideration of interaction terms in the fermionic Hamiltonian rests upon the Hartree-Fock 
procedure after fixing the appropriate gauge. Finally, we discuss the magnetic properties of this quantum 
spin model at zero as well as non-zero temperatures and analyze the validity of the approximation used 
through a comparison with the results of the exact diagonalization method for finite (up to 36 spins) 
chains. Besides the m = 1/3 plateau the most prominent feature of the magnetization curve is a jump at 
intermediate field present for certain values of the frustrating vertical bond. 
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! 1 Introduction 

■ Low-dimensional quantum spin models on frustrated lat- 
, tices represent objects of intense current research (see e.g. 
' Refs.[l][2l[3l|4] for the recent reviews). The analytical study 
\ of such models is however quite involved, since the inter- 
I play of competing interactions, quantum fluctuations and 
' magnetic field may produce at sufficiently low tempera- 
, tures a diversity of unusual quantum phases. A dimeriza- 
' tion in the ground state, localized excitations, magnctiza- 
\ tion jumps and plateaus are the most typical phenomena, 
. which might possibly appear in geometrically frustrated 
' quantum spin models[l][2l[3j|4l[5j|6]. 

, The spin-^ Heisenberg model on the diamond chain 
is an example of a frustrated spin system, which repre- 
sents one of the simplest quantum models with the exactly 
known monomer-dimer ground state. When the nearest- 
neighbor spins from a diamond chain are coupled through 
antiferromagnetic interactions, this model actually exhibits 
either ferrimagnetically ordered ground states or the disor- 
dered tetramer-dimer and dimer-monomer spin-fluid phases 
depending on the relative strength of the geometric frus- 
tration [7] . A number of exact results for low-temperature 
properties are also available at sufficiently high magnetic 
fields close to a saturation value at which this quantum 
system shows a magnetization jump towards the com- 



pletely polarized phase (see e.g. Refs. [HIS] and [S]). Exact 
results can be also found for the special Ising-Heisenberg 
diamond chain when the only quantum interaction is be- 
tween spins on vertical bonds and one may apply the 
decoration-iteration transformation to study this simpli- 
fied model rigorously in the whole temperature range [5]. 
Most recently the ground states of the mixed diamond 
chain with higher spins have been found rigorously [H]. 
Quantum phase transitions and finite-temperature prop- 
erties of the diamond chain with the mixed spins 1 and 
1/2 have been studied in Refs. jlOlfTT] . 

The ground-state phase diagram of the spin-^ Heisen- 
berg model on the distorted diamond chain cannot be 
found exactly in general, but the extensive study of the 
ground phase has been performed in Refs. P^IT^ by means 
of the exact diagonalization and some perturbative ap- 
proaches. The effect of the exchange anisotropy on the 
ground state properties of the spin-i XXZ diamond chain 
was considered in Refs. [HlfTSlfTB] . where the interesting 
inversion phenomenon has been theoretically predicted. 
In the case of the easy-plane anisotropy, the ground-state 
phase diagram contains an additional Neel phase as a re- 
sult of the interplay between the exchange anisotropy, ge- 
ometric frustration and quantum fluctuations. 
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The experimental detection of the frustrated diamond 
chain in azurite Cu3(C03)2(OH)2 [IT], the ferrimagnetic 
diamond chain in organic-radical system |18] . as well as, 
the ferromagnetic diamond chains in polymeric coordina- 
tion compounds Bi4Cu3V20i4 [13, Cu3(Te03)2Br2 [IS] 
and Cu3(N3)6(DMF)2 [H] have stimulated a number of 
experimental I5a[23l[Ml[25l[26l[2ll and theoretical [SIMIIMI 
[501151] studies over the last few years. The study of the 
field-dependent magnetization curve in the ground state, 
dynamic properties and the determination of spin cou- 
plings are the most discussed problems nowadays p7[[251 

Most of the previous theoretical treatments are based 
on numerical techniques. Hence an analytical approach 
covering the ground state as well as thermodynamic prop- 
erties is desirable. A promising method is the Jordan- 
Wigner fcrmionization of the spin degrees of freedom. Such 
an approach has been presented in Rcfs. {^[50] . However, 
the results presented there suffer from the neglecting of 
phase factors and, therefore, they cannot be considered 
as reliable (see our discussion in Sect. [3]). Note that the 
Jordan- Wigner transformation also gives a useful repre- 
sentation for the quantum models on two-leg ladders j321 
[551151] . In the case of the railroad ladder approximative 
considerations allow to describe not only the ground state, 
but also the dynamic properties [55] . 

The goal of the present paper is to describe thermo- 
dynamic and magnetic properties of the quantum spin- 
i XX diamond chain in fcrmionic language by means 
of the Jordan- Wigner transformation and to complement 
our findings by exact diagonalization data. It should be 
stressed, however, that the diamond chain model is not 
reduced after performing the Jordan- Wigner transforma- 
tion to free fermions. The XX interaction terms are not 
mapped to the two-fermion terms but they might con- 
tain phase factors whose specific form will basically de- 
pend on a particular choice of the nonlocal Jordan- Wigner 
transformation (sec the discussion on the railroad ladder 
in Ref. [55]). On the other hand, all different fermionic 
representations are connected through appropriate gauge 
transformations. Because of the phase factors, which ef- 
fectively introduce interactions between fermions and act 
as operators, the quantum diamond chain model cannot 
be treated rigorously anymore and some further approx- 
imations are required. Note that this kind of fermionic 
representation is extremely useful for studying dynamic 
properties of quantum spin-i models (see e.g. Refs. [351 
[36]). 

The outline of this paper is as follows. In Sec. [21 the 
Jordan- Wigner transformation for the spin-i XX model 
on the diamond chain will be considered. The free-fermion 
model and its relation to the quantum spin model is an- 
alyzed in Sec. [5] The Hartree-Fock approximation is ap- 
plied to the fermionic counterpart of the symmetric and 
distorted diamond chain in Sees. H] and [SJ The ground- 
state properties will be finally calculated and compared 
with the known exact and numerical results there. The 
thermodynamic properties will be considered in Sec.[n]and 
the obtained results are summarized in Sec. [7] 




Fig. 1. Different ways of assigning Jordan- Wigner transforma- 
tion on the diamond chain are marked by thick lines. The case 
i) (panel a) and the case ii) (panel b), see the main text. 



2 Jordan-Wigner transformation for the 
spin-i XX diamond chain 

Consider the quantum spin-^ XX model on the diamond 
chain (see Fig. [1]) with the following Hamiltonian: 



N 



1=1 



+M4,i''3,i + 4r'^i,/+i) + h-c.) 



p=i 



(1) 



Here, N is the number of unit cells, s^, = , ± ist , are 



the spin raising and lowering operators, s. 



3+ S- - i 
^p,l°p.l 2' 



; (ot = X, y, z) are the usual Cartesian components of 

the Pauli spin-i operator with the first index correspond- 
ing to a sublattice and the second to a cell, h is the ex- 
ternal magnetic field (we set g^B = !)• We will further 
distinguish between two models: the distorted (Ji ^ J^) 
and the symmetric ( Ji = J3) diamond chain. 

The Jordan-Wigner transformation can be unambigu- 
ously defined on a linear chain [571155] , where all sites can 
be enumerated subsequently. The diamond chain is how- 
ever a three sublattice model and there exists at least two 
identical ways of arranging its sites in a one-dimensional 
sequence as shown on panels (a) and (b) in Fig. [T] The 
case i) sets the following order of sites . . . , {3,1 — 1), {1,1), 
(2, /), (3, /), (1, / + 1), ... and the case ii) corresponds to 
another choice of order . . . , {2,1 — 1), (1, 1), (3, 1), (2, 1), 
{1,1 + 1),.... 

Following the usual procedure one can define the non- 
local Jordan-Wigner transformation, which introduces new 
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Fermi operators by multiplying the spin lowering and rais- 
ing operators by the Jordan- Wigner factors (— 2sf^) re- 
ferring to all preceding spins {571155] . For the case i) we 
have: 



'2,/- 



=ai,;exp 



=a2,;exp 



=a3,;exp 



3 ;-i 

p—1 2—1 

3 /-I 

/ 4- . X 

1 J. ,t ^ ^ 

\ p—1 i—1 

/ 3 l-l 

i- + _i 

2 ^ '^iij^nd 2" 



(2) 



Subsequently, the spin interaction terms to emerge in the 
Hamiltonian ((T]) can also be rewritten by the use of Eqs. 
([2]) into the fermionic representation; s^jS^j = a^ja2,/, 

-t- — + -t- — + i-nat ,a-> i 

+ — + ITT at , a I -L, _ _L 

^2,1^1,1+1 = Oj^^e 3,, 3''ai,/+i, S3^Si,_^i = a^iUi^i+i. 
Hence, it follows that the fermionic representation of the 
Hamiltonian of the spin-^ XX diamond chain then reads: 



1 ^ r 

Hxx = 2^ {J2a'2.i^i.i + Ji{0'i^iO'2d + a^^iaijj^i) 
1=1 L 



-h.c.) - 2/i^ a+iGpj - - 



p=i 



(3) 



In the current ordering the sites (1,/), (3,Z) and (2,1), 
{1,1 + 1) arc not nearest neighbors. Therefore, the cor- 
responding couplings contain phase factors in the hop- 
ping terms that prevent a rigorous consideration of the 
investigated model. By imposing periodic boundary con- 
ditions in the considered model isf^ transform into 
many fermion terms. Such boundary terms may be ne- 
glected after making the Jordan- Wigner transformation, 
since they are irrelevant for the study of static properties 
in the thermodynamic limit |37ll38j . 

Considering the numbering of sites ii) the Jordan- 
Wigner transformation should be defined as: 



=ai ;exp 



3 l-l 



0=1 i=l 



S2 z=a2,;exp 



'3,r 



=03 ;exp 



(3 l-l 
p=l i=l 



3 l-l \" 
p=l i=l / . 



(4) 



Hence, one gets another fermionic representation of the 
original Hamiltonian ([T]): 



1 ^ 



1=1 



{,720.2 ids^i + J3(a^ ,03^; + ;ai^;+i) 



+Ji(,a^;e ' a2,/ + flg ;e ' ai_/_|_ij 



-h.c.) -2/7^ (d+idp,i - i 
p=i ^ 



(5) 



It should be noted that the Hamiltonian ([5]) differs from 
the Hamiltonian (|3]). The hopping terms, which have not 
contained the fermion interaction before, now acquire a 
phase factor, and conversely, in other terms the relevant 
phase factors have disappeared. A common feature of both 
Hamiltonians is that they contain the interaction only at 
some particular bonds, and therefore, if one intends to 
treat the model approximatively, the contributions of dif- 
ferent bonds will not be generally considered on equal foot- 
ing. 

Using (0), (H]) one can check that operators ap^m and 
o-p,m are connected through the gauge transformation: 

ai, 

a2d=d2,i exp {iTTa^iUs^i), 0,2^1=02,1 exp {iira^ 103,1), 
03,1=03,1 exp {iTTO^ [02, i), 03,1=03,1 exp {iTra^^a2,i)- (6) 

Thus, the Hamiltonians ([5]) and ([5|) describe the same 
physics. 

One can also consider a more general gauge transfor- 
mation for Fermi operators Cx- 



Cx = c'^y*^ "'"y^-^Cx, 



(7) 



where x, y are general indices which may include the site 
and sublattice numbers. If axy = ayx + ^im new operators 
Cx satisfy the Fermi commutation relations. 

Fixing the gauge in transformation ([7]), we can con- 
sider a particular case: 



Cl,l = Oi,l, 



C2,i = e 



'02,;, 02,; = e '^"^■'"''''€2,1, 



C3^i = e^2-.,:-2.^a3,i,a3,i = e-^^^=>'^=''c3,,. (8) 

This gives as a result the Hamiltonian in the most sym- 
metric form: 



Hx 



(Ji(c+ e 



C2,l + ;e 



1 ^ 

-y 

2 ^ 

1=1 



-HJ3(c+e^^=2.'^^-'c3,/ + c+e-'^=3,,-3, 



'ci,i+i) 



Cl,l+lj 



+ J2C2,lC3,l + ll-C.) - 2/1^ (cp^lCp,l - ^ 



p=l 



(9) 
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One can also easily find the spin-fermion transforma- 
tion that leads to the representation It takes the form: 



not invariant under such transformation because all the 
terms change their signs: 



^2 ;=C2,/exp 



3 i-1 
p— 1 i—1 



N 



(10) 



3 /-I \ ■ 

P=l 4=1 / . 



^3,r 



=C3,zexp 



V p=i i=i 



Finally, one should note that equations (|10p represent a 
more general formulation of the Jordan- Wigner transfor- 
mation for the spin ladders: 



(11) 



with the condition ipx.y = ify.x^^ [33j . Here x, y are again 
general indices which may include the site and sublattice 
numbers. 

Wc would like also to comment some limiting cases 
of the model. If Ji = or J3 = the model resembles 
the exactly solvable trimerized XX chain [551l40lHT| . In- 
deed, if J3 = (or Ji = 0) wc come after the Jordan- 
Wigner transformation ([2]) (or (|4])) to the Hamiltonian 
^ (or ([5])) in the form of free fermions. In case J2 = 
but Ji 7^ 0, J3 7^ it is not possible to find such a form of 
the Jordan- Wigner transformation that avoids the inter- 
action between fermions in the transformed Hamiltonian. 
Therefore, the consideration of the models with Ji 7^ 0, 
J3 7^ requires further approximations even if J2 = 0. 

For the further consideration it is useful to remind that 
the I Dn (I t)n) spin state corresponds to the empty |0)„ 
(filled |1)„ = c+|0)„) fermion state. Thus, the magneti- 
zation per spin in z-direction is related to the fermionic 
averages in the following way: 



^ 3 N ^ 3 N 



3N 



p=i 1=1 



=1 ;=i 



(12) 



" 1=1 

+J3{at,i'^3,i + aj/ai,i+i) + h.c.^ . (13) 

One may also use a kind of phase transformation (ap_/ = 
(— l)ap_/) for some sites to revert the signs. However, five 
terms of the Hamiltonian p3D should reverse their signs 
for every particular I using transformations for only three 
operators and hence it follows that it is impossible to 
satisfy simultaneously all the conditions quite generally. 
However, such a possibility appears for some particular 
cases Ji = or J3 = which correspond to the exactly 
solvable trimerized chain [39ll40ll4T| . The situation for the 
Hamiltonians ([3]) and ([SJ is different. After the particle- 
hole transformation the terms which contain phase factors 
preserve their signs, the signs of other two-fermionic terms 
can be always reverted by the aforementioned phase trans- 
formations. 

Thus, the symmetry of the spin-i model is broken in 
the free-fermion Hamiltonian. It may mean that the mag- 
netization calculated from this Hamiltonian is not an odd 
function of the external magnetic field anymore. We can 
prove it in a way similar to [42] . One has to calculate the 
magnetization using the fermionic representation: 

where /3 — 1/T is the inverse temperature. We have also 
introduced the density of states 



3 N 



(15) 



Ap_i denote the eigenvalues of the free-fermion Hamilto- 
nian (|13|) . Generally, one has to prove the asymmetry of 
the density of states for h = Q. The simplest way is to 
show that the odd moments for the density of states 



3 Phase factors and their importance. 
Free-fermion model on diamond chain 

Now we would like to analyze the general properties of the 
free-fermion model on the diamond chain, i.e., the model 
which one gets after the Jordan- Wigner transformation 
ignoring the phase factors terms [53] . 

The Hamiltonian of the spin-i XX model, as well as a 
more general XXZ model, in zero field is invariant under 
the reflection of z-componcnt of all spins: s 



<i [I t 

/p,i ^ I fermionic language it corresponds to the 

particle-hole transformation: Up^i -> a^j [|0)p,; — > |l)p,;]. 
It is clear that the Hamiltonian of the free-fermion gas is 



M 



dEE^p{E) 



can take non-zero values. Using the Green's function method 
1431. we find the third moment in the form: 



1 



h{Jl+2J'l+2j'. 



1 



JlJ2J3- (16) 



Hence, it remains nonzero in the zero-field case and gives 
a nonzero contribution to the zero-field magnetization. 

This general conclusion can be confirmed by numer- 
ical calculations. The free-fermion Hamiltonian (jl3p can 
be diagonalized by subsequent application of the Fourier 
transformation (which present the Hamiltonian as a 3 x 3 
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Fig. 2. (Color onhne) The ground-state magnetization vs. the 
external field of the free-fermion (FF) model on: a) the sym- 
metric diamond chain Ji = Js = 1, J2 = 0, 1, 2, the thick 
step-like lines correspond to the exact diagonalization (ED) of 
the spin model ([T| of 36 spins for J2 = 0, 1; b) the distorted 
diamond chain J2 — 2, Ji = 1, J3 = 0, 0.5, 1, the thick step-like 
line corresponds to the exact diagonalization of the spin model 
of 36 spins (m for J3 = 0.5. 



quadratic form) and unitary transformation. The eigenen- 
ergies ^„(k) = h + w„(k) are expressed through the solu- 
tions lo„{k) of the cubic equation: 



l^n(K)^ +p{n)^^n{K) + !?(«) = 0, 



p{k) 



Jl Jl Jl 



J2 



J1J3 COS K, 



qin) = -^(2JiJ3 + {Jf + Ji)cosK.). 

Here k is the wave vector of the Fourier transformation, 
and n is the sublattice index. The density of states is 
rewritten as p{E) = J^n Pn{E) where Pn{E) = Y,,, ^{E - 
h — ll)„(k)). Using this representation of the density of 
states and the expression for the magnetization ([H)) one 
can easily understand the curves shown in Fig. [21 The gap 
between nearest bands corresponds to the plateau of the 
magnetization. The jumps of the magnetization is a mark 
of the flat band in the model. Indeed, for the symmetric 
chain one can find the solution in simple form: uji{k) = 



-J2/2, L^2,3(«) = J2(l ± ^/l+{Ui/J2)Hl + COaK))/4. 

Thus, the ground-state magnetization has the jump at 
h = — J2/2. Note that in case all couplings Jp are non- 
zero, the magnetization gains a finite value even for zero 
fields (see Fig. , which is clearly an artifact of the used 
free-fermion approximation. In Fig.[5]we also compare the 
magnetization derived from the free-fermion model (1131) 
with the exact diagonalization data of model (H} . One can 
see that the good agreement is achieved for the strong neg- 
ative fields only where the system is close to the ordered 
state with all spins down. 

To conclude, the free-fermion model on the diamond 
chain generally looses the symmetry of the corresponding 
quantum spin model and, moreover, it gives physically not 
acceptable solution for the nonzero magnetization even in 
the zero-field case. 



4 Hartree-Fock approximation for the 
distorted diamond chain 

When the sites (2,1) and (3,1) are shifted from its sym- 
metric positions, we come to a distorted diamond chain 
where only some of the spin interactions differ from each 
other. Without loss of the generality one may impose that 
Jl > J3 and choose such a fermionic representation of the 
Hamiltonian ^ that the fermion interaction terms will 
finally appear along the weaker bonds only. The latter 
terms can be then treated perturbatively. Using the alge- 
bra of Fermi operators, i.e., exp{inap [Up^i) = 1 — 2a^iap^i, 
one can rewrite ^ in the following form: 



1 ^ 

2 ^ 



+ J3(a+,(l-2aJ;a2,/)a3,/+a+,(l-2a+ja3,;)ai.;+i) 

^ / 1 " 

+ h.c.)-2h'^{a+iapj-- 



(17) 



Applying the Hartree-Fock or mean- field- like approxima- 
tion which preserves the correlations between 
neighboring sites, e.g. 



-a+;a3,/Ai - ^30+, 02,; + A3A1 
+a^ja2jn2 + A2a^ja3,i ~ ^2n^2, (18) 



Hx 



we come to the quadratic form in Fermi operators: 
1 ^ 

2E[((^2+4J3A2)a+,a3,/ 

-|-(Ji + 2J3Ai)(a+,a2,; -f a+jai,;+i) 
-1-/3(1 - 2n2)(a;|''_;a3^; -f aj;ai,/+i) + h.c.) 

-2ha\^ax^i - 2(/i + .hM){a^^a2,i + a;^/a3,/) 
+iVeo. (19) 
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Here we use the following notations: Ai = {a^ 10.3,1), A2 = 

n2 = (aj;a2,;) = {at,ia3,i), eo = ^ - 4J3A1 A2 +4J3n2A3. 
In the elementary contractions A2, A3 the invariance of 
the initial Hamiltonian with respect to the space reflec- 
tion was exploited. 

After Fourier transformation, 

^ — 1 

iV iV 
K = 2'Km/N, m = — — + 1, . . . , "^(for even A^), 

iV - 1 iV - 1 
m = ^ — , . . . , — - — (for odd A), 

the Hamiltonian can be represented in the matrix form 
3 

Hxx=Y. '^pA^H.^ag,, + Neo, (20) 

K, p,q—l 

Hn{K) = -/i,%2(k) =H33(k) = -h-2A3J3, 
Hi2(«;)=H;i(K)-i(Ji+2J3Ai+J3(l-2n2)e-'''), 

Hi3{n)=Hl,{K)^^iJ3{l-2n2) + {Ji+'^J3A,)c-n, 

H23(Ai) = n32iK) = i(J2 + 4^2 J3)- 

This Hamiltonian can be reduced to a diagonal form: 

H^^ = ^ Ap{K)Tj+^T]p,^ + Neo, (21) 

by some unitary transformation rip,^ ~ X]g=i '^p.gC'^)'^?-^ 
where Up^qin) is a unitary matrix. Generally, it corre- 
sponds to the 3x3 eigenvalue and eigenvector problems. 
Finally, we approximately presented the model as the free- 
fermion gas, and its thermodynamic and correlation func- 
tions can be calculated straightforwardly. However, the el- 
ementary contractions, which are included into the Hamil- 
tonian parameters, are unknown and have to be found 
self-consistcntly: 

n2==(a+;a2,/)=^^|u2,p(K)P(?y,!;K'7p,A.), 

Ai = (a+ 03,/)= ^ ^ ^,p('«)"3,p(k)(?7^«'?p,«), 

K, p 

A3^{a+ia3j)=J2'Yui,pi^)u3A'^){Vp,KVp.K)- (22) 

Here (...) means the thermodynamic average with the 
Hamiltonian (PT|) . and for the ideal Fermi gas (??jJ'k^p,k) = 
l/{exp{(3Ap{K)) + 1) is the Fermi function. Thus, the ini- 
tial statistical mechanical problem for the spin model is 




h 



Fig. 3. (Color online) The ground-state magnetization vs. 
the external field of the distorted diamond chain for Ji — 1, 
J3 = 0,0.25,0.5,0.75,1. a) J2 = 1. Inset shows the com- 
parison of the results of the Hartree-Fock approximation and 
exact diagonalization for A'^ = 24, 30, 36 (step-like lines) for 
J3 = 0.25, 0.5. b) J2 — 2. Inset shows the comparison of the 
results of the Hartree-Fock approximation and exact diagonal- 
ization for A'^ — 24, 30, 36 (step-like lines) for J3 = 0.5. 



turned into the solution of the set of equations (P^ . We 
have solved this problem numerically and then used the 
formulae for the free-fermion gas to calculate the thermo- 
dynamics and static properties of the model. Particularly, 
the equation (fT4)) can be used for calculating the total 
magnetization. 

For T = {Vp.kVp,^) becomes the Heaviside step func- 
tion and the summation in (j22p can be restricted to us 
which satisfy Ap{K) < 0. This simplifies the calculation 
of the ground state properties of the model. The results 
for the ground state magnetization are shown in Fig. [3] 
It is worthwhile to remark that one recovers the exact 
result for the magnetization curve of the uniform j371l38] 
and trimerizcd [iOlHT] spin-i XX chains in the partic- 
ular cases Ji = J2 = 1, J3 = and Ji = 1, J2 = 2, 
J3 = shown in Fig. [3] Another particular case, which is 
amenable for an exact calculation, is the dimer-monomer 
limit J2 = 2, Ji = J3 = 1 (see Fig.[3jD). The magnetization 
has a step-like shape in this special case (see Appendix). 
The Hartree-Fock equations have the following solu- 
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tions for the elementary contractions: 712 = ^, = — ^, 
^2 = ^3 = for ft, < Inserting the solutions obtained 
for J2 = 2, Ji = J3 = 1 into Eq. (fTO)) . one comes to the 
approximative Hamiltonian in the form 



1 ^ 



1=1 



-2h [ c+ ci,; + c+;C2,; + c+;C3,; - - 



(23) 



It can be diagonalized by the canonical transformation: 

Vs,l = ^(C2,i-C3,i), r]t,l = -^{C2,l+C3,l), VmJ = cij. As a 

result we come to the free-fermion gas of the form: 



1 ^ 



'^hinliTlt.i+vXiVs.i+Vm.^iVm.v 



(24) 



The ground state of the model can be easily found: \GS) ~ 



bands separated by gaps. The external field plays the role 
of the chemical potential of fermions and thus, it can shift 
their Fermi level to the area of a gap between two bands. 
Such a situation corresponds to the 1/3-plateau of the 
magnetization. It can be also understood that all elemen- 
tary contractions do not change along the plateau. There- 
fore, the width of the plateau is equal to the width of the 
gap in case the magnetization curve does not reveal any 
jumps. 

The effect of J3 interaction on the ground-state mag- 
netization process can be understood from Fig. |H We see 



(a) 



riz 7i(cL - c+ )io) = ri/i i)iAWA if ft < 0, 

\GS) = n,7lC+«, -c+)|0) = a I t)i,a2,/;3,?] if 
ft > 0. Here [2,Z;3,Z] represents the singlet dimer state 
between spins on sites (2,/) and (3, Z). As one can see, 
the Hamiltonian ((24|) indeed represents the model of non- 
correlated singlet dimers created between the spins on 
sites {2,1) and (3,Z), which are separated through the 
free monomeric spins residing the sites (1,0 snd thus, 
there cannot be any correlation between singlet dimers. 
We note that the same approach also recovers the com- 
pletely dimerized state in the Majumdar- Ghosh limit of 
the zig-zag ladder [51] . 

For the intermediate values of J3 only approximative 
results are available. To understand to what extent the 
elaborated approach is valid we compare it with the re- 
sults of the exact diagonalization of finite chain up to 36 
spins. For some particular values of J3 such comparison 
is shown in the insets to Fig. |21 It can be seen that the 
approximative results are quite accurate for sufficiently 
strong fields ft, and one recovers the 1/3-plateau for the di- 
amond chains. Therefore, the Hartree-Fock method is reli- 
able even for the case of strong frustration ( J3 ~ Ji). How- 
ever, for small field it may produce an artificial jump of the 
magnetization for J3 = 0.25, 0.5, 0.75 if one compares the 
results stemming from the Hartree-Fock approach with 
the data of the exact numerical diagonalization method 
(see. Fig.|3ji). 

The comparison with the free-fermion model, see Sect.|31 
demonstrates that neglecting of the phase factors leads to 
quite different magnetization curves. Besides the wrong 
low-field behavior there is also a significant difference in 
the saturation field. 

The 1/3-plateau in the magnetization curve can be 
understood from formula (|14p as an integral over the den- 
sity of states. The approximative representation of the ini- 
tial spin model ([l} is a fermionic model with three types 
of fermions. Thus, the energy spectrum consists of three 
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Fig. 4. The ground-state phase diagram of the distorted di- 
amond chain in the Hartree-Fock approximation. The upper 
(lower) panel corresponds to the case Ji = 1, J2 = 1 (Ji = 1, 
J2 = 2). 



that for initially uniform chain (J2 = ^/i = 1) an addition 
of the frustrating interaction J3 generally causes an onset 
of the intermediate magnetization plateau at m = 1/3, 
which becomes the broader the stronger the interaction 
parameter J3 is (see Fig.HJi). Generally, the coupling J3 in- 
creases the width of the magnetization plateau also in case 
of the initial trimerizcd chain. It is interesting to note that 
the upper critical field for the 1/3-magnetization plateau 
in case Ji = 1, J2 = 2 does not show any marked depen- 
dence on J3. 
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It should be noted that the ground state phase diagram 
of the generally anisotropic distorted diamond chain was 
studied previously P^ITSIITC] , The so-called inversion ar- 
eas may be found in this more general model where the 
Neel phase becomes the ground state in the parameter 
space with the predominant XF-like interaction and the 
spin-fluid phase contrary appears in the parameter region 
with the predominant Ising-like interaction. The Neel or- 
der considered in Refs. [Hl[T^[TB] is related to the antifer- 
romagnetic ordering of spin blocks (s3.;_i -I- Si^; -I- 82,;), 
and thus corresponds to 6 spins in the magnetic cell. This 
phase seems to disappear in the pure XY limit (see Fig. 2 
of Ref.[l6]V at least it is true for Ji S> J2, J3. Since we 
have not considered the possibility of doubling of the el- 
ementary magnetic cell in our approach, we cannot ver- 
ify whether such a phase exists. Moreover, the spin-fluid 
phase found in Ref. (TH] for ft, = we do not find in our ap- 
proach, which might be attributed to its mean-field char- 
acter. 



5 Hartree-Fock approximation for the 
symmetric diamond chain 

Here we consider the case of symmetric diamond chain 
when the interactions Ji = J3. For this purpose, let us 
consider the fermionic representation Q of the spin Hamil- 
tonian which provides a symmetric representation for those 
bonds. In the fermionic Hamiltonian we keep all the terms 
that correspond to the correlation between neighboring 
sites in the decoupling. Since the coefficients of the Hamil- 
tonian ^ are complex the elementary contractions can 
also be complex valued. Similarly to previous section we 
introduce the following notations: 

A3 - (c+,C3,z) = (c+ +iC2j). (25) 

Here we additionally used the invariance of the initial spin 
Hamiltonian with the respect to the space reflection. The 
invariance to the exchange of the second and third sublat- 
tices, gives us the following connections between different 
contractions: A2 = A*^ {A^ = A§, A{ = -A^). Here 
(j4p) is the real (imaginary) part of Ap. The Hamiltonian 
within the adopted approximation will be as follows: 

1 " 

H.x = - Y.'^J2 + 4Ji(Af -f ^^)](C+ C3,, + C+ C2,0 

1=1 

+Ji («( + <i+i)(c2,/ + C3,0 + h.c.) 
+iJi {{ct^i+i - c^;,)(c2,/ - C3,/) - h.c.) 
-2(/ic+ ci,,+/i2C+ C2j+/i24,C3,i)] +iVeo (26) 

where ^ Ji(l - 7I2 + ), ]{ = Ji(n2 + ), ^2 = 
h-2MA^-Ai), eo = f -4JiAf (A«+A^)+4Jin2(Af- 
Ai). 



After performing Fourier transformation, the Hamilto- 
nian can be written in a matrix form as follows 

3 

Hxx = X! X! ^P.9('*)"^«'^9.'= + ^'^0 (27) 

K. p.q — 1 

Uii{k) = -/i, 

H22('t) = H33('*) = -ft - 2Ji(A^ - Ai), 

^i3(^)-H5i(^)-iv'('^i'')'+(^~i')'(°'^+^"'^*^'^^)' 

H23(a*) = H32(a«) = \{J2 + 4Jl(A« + Ai)), (28) 

where Ji = y^(Jf)2 -|- (j/)2, tan0 = J^/Jf . 

Similarly to the previous section the problem is rewrit- 
ten as the free-fermion gas. The Hamiltonian (j26p or ([^7)) 
contains the unknown contractions ri2, Ai, A2. These con- 
tractions have to be found from the set of self-consistent 
equations identical to ((22|) . 

The results for the ground-state magnetization are pre- 
sented in Fig. O We set Ji = 1 here. One can notice that 
the magnetization for J2 = 2Ji corresponds to the ex- 
act one (see Appendix). Following the arguments of the 
previous section, it can be shown that the Hartree-Fock 
approximation leads to the dimer-monomer ground state 
which is exact for this case. For J2 > 2Ji the magneti- 
zation curve still has a step-like shape and corresponds 
to the exact result. For intermediate values of J2 we also 
obtain rather good coincidence with the exact diagonal- 
ization data. We even recover the magnetization jump for 
intermediate fields which is present also in the exact diag- 
onalization data. However the approximate results for the 
low-field magnetization show a noticeable difference with 
the exact diagonalization data. Moreover for J2 > O.86J1 
we get an incorrect non-zero magnetization even for zero 
fields ft. It is well known that if 0.909 < J2/J1 < 2 the 
singlet dimer-tetramer phase is the ground state of the 
isotropic diamond chain [7] . This phase survives also in the 
model with the easy-plane anisotropy |14j . The elemen- 
tary magnetic cell is doubled in this case. Therefore, the 
current approach, where a doubling of the unit cell is not 
taken into account, does not capture this phase. At higher 
fields the dimer-tetramer order is destroyed and the order- 
ing of the singlet dimers on vertical bonds becomes more 
favorable (see Fig.[S}D). Since the elementary magnetic cell 
is not doubled anymore, we obtain an excellent agreement 
between the approximate and exact diagonalization data 
(see the inset in Fig. [5^) . Our exact diagonalization data 
are also in agreement with the previous results for the 
easy-plane XX Z diamond chain P^[T51[T5] . Note that the 
model with the easy-plane anisotropy in zero field is char- 
acterized by the spin-fluid phase for J2 < Ji (see Fig. 3 in 
Ref. [13]) in contrast to the isotropic model where the ferri- 
magnetic order persists in the ground state [7] . The exact 
diagonalization data show the continuous growth of the 
ground-state magnetization with field for J2 ^ Ji- This 
is usually the sign of a gapless excitation spectrum and 
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(a) 0.5 



(b) 0.5 




Fig. 5. (Color online) The ground-state magnetization vs. 
the external field of the symmetric diamond chain: Ji — 1 
a) the results for the Hartree-Fock approximation are shown 
and compared with the exact diagonalization data in the inset 
for J2 — 0, 1, 1.75; b) the results of the exact diagonalization 
for finite chains of A'^ = 30, 36 spins with J2 = 0, 0.5, 1, 1.75. 



the spin-fluid phase. For stronger J2 couplings we obtain 
the zero-magnetization plateau and a very steep increase 
of the magnetization to the 1/3 plateau for some criti- 
cal field {h sa O.I68J1 for J2 = 1.75Ji). The phase with 
the zero magnetization plateau can be identified as the 
singlet dimer-tetramer phase, and critical field as a field 
which destroys tetramers in the ground state. Consider- 
ing the ground-state phase diagram presented in Fig. |6l 
we have to conclude that the considered Hartree-Fock ap- 
proach is not able to exhibit the singlet dimer-tetramer 
phase with the zero magnetization plateau which appears 
in the thin region of small fields and 1 < J2/J1 < 2 in 
accordance with the exact diagonalization data. However, 
as soon as the magnetic cell is not doubled anymore, the 
ground state phase diagram are in good agreement with 
the exact diagonalization data. 

Again, the comparison with the free-fermion model, see 
Sect. [21 shows that it is not a reasonable approximation 
to neglect the phase factors. In particular, the jump found 
within our approach and confirmed by exact diagonaliza- 
tion is not present in the free-fermion model. Comparing 
the results of the previous section we see that both for- 
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Fig. 6. The ground-state phase diagram in the h — J2 plane 
for the symmetric diamond chain with Ji = J3 = 1 in the 
Hartree-Fock approximation. 



mulations with different gauges produces the same results 
for the ground state properties. Thus the Hartree-Fock 
approximation preserves the gauge invariance. 

The results of the exact diagonalization show also some 
indication of the second plateau at 2/3 of the saturated 
magnetization. We note that the elaborated Hartree-Fock 
approach cannot reproduce such a behavior. To get it one 
has to consider the possibility of doubling of the elemen- 
tary magnetic cell in the decoupling procedure which was 
not assumed in Eqs. (^5)) and 



6 Diamond chain at non-zero temperatures 

In this section we consider the thermodynamic proper- 
ties of the diamond chain within the fermionization ap- 
proach and the Hartree-Fock approximation. After hav- 
ing rewritten the spin Hamiltonian as a quadratic form of 
Fermi operators, one can find its eigenvalues and obtain 
the Helmholtz free energy per site: 

p=l K 

1 r°° 

= — dEpiE)He-P'' + l), (29) 

P J — 00 

the internal energy per site 

u=^{H)^ I dEp{E)n{E)E + eo, (30) 

and the specific heat per site c = ^ 

run 2 

du dn2 



rt— 1 ^ 

— run ^ 



Here /? = 1 /T is the inverse temperature. 



dT ■ 
(31) 
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Fig. 7. (Color online) The specific heat in the Hartree-Fock ap- Fig. 8. (Color online) The specific heat in the Hartree-Fock 
proximation vs. temperature of the symmetric diamond chain approximation vs. temperature of the distorted diamond chain 



(Ji = Ja = 1), Ja = 0,0.5, 1.5,2: a) h = 0, b) /i = 0.1. 



Ji = 1, J2 = 2, Js = 0, 0.25, 0.5, 0.75, 1: a.) h = 0, h) h = 0.2. 



To calculate the specific heat we perforin the numerical 
difFerentiation. The results of computations are shown in 
Figs.[7]and|51 Some typical features can be observed there. 
The specific heat for small vertical couplings J2 grows up 
linearly with the temperature which is inherent for sys- 
tems with a gapless excitation spectrum. For sufficiently 
strong J2 (see the curves corresponding to J2 = 1.5,2 
in Fig. [T]) the ground state is degenerate due to dimer- 
monomer phase, and the system has a gap between the 
ground and excited states. Therefore, we observe the ex- 
ponential growth of the specific heat with temperature. 
When a small external field h is applied, it does not pro- 
duce a gap in the excitation spectrum, thus the tempera- 
ture dependence of the specific heat remains qualitatively 
the same. The effect of the external field on the dimer- 
monomer phase is somewhat different. The ground state 
degeneracy is lifted due to the Zceman term —hj^p i Sp 
As a results we have two kinds of excitations with differ- 
ent gaps. This leads to the appearance of the double-peak 
structure in the temperature dependence of the specific 
heat. The first peak corresponds to the low-energy flip ex- 
citations of the monomer spins si^i, and the second one 
corresponds to the thermal excitations of the spins cre- 
ating singlet dimcrs. Here we should emphasize that the 



dimer-monomer phase as well as the ground state degen- 
eracy is an artifact of the elaborated approximation for 
0.909 < J2/J1 < 2 (see discussion in Sec. 5). Therefore, 
the results presented in Fig. [7] for J2 = 1-5 cannot repro- 
duce the features of the specific heat of the spin-i XX 
diamond chain. 

In contrast to the symmetric chain, the distorted dia- 
mond chain shows the double-peak structure in the tem- 
perature dependence of the specific heat even in the zero- 
field case (see Fig. , whereas the external field can even 
demolish this double-peak structure. The mentioned fea- 
tures are not only the consequence of the XX anisotropy 
in the spin interaction. The similar behavior of the temper- 
ature dependent specific heat with the double-peak struc- 
ture was also reported for the Heisenberg diamond chain, 
see Fig. 8 of Ref. ^U\. The external field can also lead 
to less visible double-peak structure in the specific heat of 
the Heisenberg diamond chain [15] . 

The field dependent magnetization for non-zero tem- 
peratures can be also understood. The plateaux and jumps 
of the magnetization are smeared out by temperature (see 
Fig.[9l). One exception is the model when J2 ^ Ji when the 
mean-field treatment preserves jumps even for small tem- 
peratures. Thus, the Hartree-Fock approximation turns 
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Fig. 9. (Color online) The magnetization in the Hartree-Fock 
approximation vs. the external field of the diamond chain for 
T = 0.05 Ji: a) Ji = Jg = 1, J2 = 0, 0.5, 1, 1.75, 2; b) Ji = J2 = 
1, Js = 0,0.25,0.5,0.75, 1 

out to be incorrect for Ji ^ J2 ^ J3 and small exter- 
nal field h, since the spontaneous (i.e. non-zero) mag- 
netization in the zero-field limit is prohibited for one- 
dimensional system at any finite temperature due to the 
Mermin- Wagner theorem [49] , 



7 Conclusions 

In the present work, we consider the Jordan- Wigner trans- 
formation for the spin-i XX model on the diamond chain 
by assuming a quite general distorted case. 

At first we showed that the free-fermion model on the 
diamond chain cannot describe the spin system since it 
loses the symmetries of a spin model and exhibits a non- 
zero magnetization in zero magnetic field. In the case of 
the distorted chain we apply the Hartree-Fock approxima- 
tion for the fermionic representation of the model which 
considers the fermion interaction along weaker bonds. For 
the symmetric diamond chain we suggest the generaliza- 
tion of the Jordan- Wigner transformation and build a fully 
symmetric fermionic representation of the spin- model. Fur- 
ther, we use the Hartree-Fock approximation to study the 
model near the limit of the dimer-monomer phase when 



the correlation between the spins from the vertical bonds 
are the strongest one. The Hartree-Fock approximation 
reproduces an exact result for the dimer-monomer ground 
state of the symmetric diamond chain (J2 > 2Ji). More- 
over, we have found that the solutions of the Hartree- 
Fock approximation for the symmetric diamond chain are 
invariant with respect to the gauge transformations. Our 
results show also good agreement with the exact diagonal- 
ization data reproducing the magnetic properties at high 
fields or small frustrations. Summarizing our findings, the 
elaborated approach for the XX diamond chain repro- 
duces a plateau in the magnetization curve at 1/3 of the 
saturation magnetization and an additional peak in the 
specific heat curve. In fermionic language a 1/3 plateau 
is caused by the gap between two fermionic bands. When 
the gap becomes large enough, the temperature depen- 
dence of the specific heat gains a distinct two-peak struc- 
ture. Note that both features are typical for azurite P51l25l 
[^[?7] although the Heisenberg model seems to be more 
appropriate for that compound. 

We have also observed some drawbacks of the Hartree- 
Fock treatment. It becomes invalid for Ji ^ J2 ^ J^- In 
this case mean-fields leads to the non-zero magnetization 
in zero fields, and the magnetization jump survives also 
for small temperatures. 

The authors would hke to thank O. Derzhko for tire useful 
discussions. T.V. was supported by the National Scholarship 
Programme of the Slovak Republic. J.S. and M.J. acknowledge 
the financial support by the ERDF EU (European Union Eu- 
ropean regional development fund) grant provided under the 
contract No. ITMS26220120005 (activity 3.2.). 

Appendix: Exact monomer-dimer ground state 
for the anisotropic diamond chain 

In case Ji = J3 > and J2 > 2 Ji the ground state of the 
Hamiltonian on the diamond chain is the dimer-monomer 
state [7]. To prove this we can follow the arguments of 
Shastry and Sutherland to obtain the ground state [TT] , 
We have to consider the Hamiltonian as a sum of Hamil- 
tonians of triangles: 

I 

Hl.l = ■yl^a.i • S3,i)zi -I- Ji ((Si^; • S,2^i)a 

+ (si,Z-S3,/)^)-^(si,/ + S2,z + S3,i)' 

-f (Si,z+i • S3,z)zv) - {sl.i+i + si I + sl i) . 

Here we introduced the notation (sp_; • Sq_m)A = Sp i^q^m + 
Sp jSq^m + i^q,m- The parts of the systems described 
by Hij and Hi.i+i are topologically identical. Therefore, 
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they have the same eigenvalues and analogous eigenstates. 
One can use the variational principle that implies Eq > 
N{ei + ei), where e; is the ground state energy of Hi i or 

Let us consider the case J2 > 2Ji. The direct calcu- 
lation shows that if < /i^ = ^(J^+^^^^^+J^ the dimer- 
monomer state |ai,/)[2^, 3/] of Hi i or is the state 

with the lowest eigenvalue ef™ = _(A+^}J2_ _ 
Sfgf. = ±i. Here \ai^i) denotes the spin up state | fi,;) 
if /i > 0, or the spin down state | lij) if ft, < 0. Thus us- 
ing the variational principle the ground state of the whole 
crystal can be built by translation of the state \aij)[2l, 31] 
over all sites. The total magnetization is ^sgn(/i) and the 
magnetization per spin is 1/6. Contrary, ii \h\ > he, the 
completely polarized state with Sf^^ = if becomes the 
state with the lowest eigenvalue e^^^ — ■^(■^^+"^•^1) _ 
Thus, the total ground state is the state with all spins up 
or down, with the magnetization per spin = 1/2. 

To summarize, the field dependence of the ground- 
state magnetization has a step-like form for J2 > 2 Ji, i.e., 
771^ = sgn(/i)/6 if \h\ < he and = sign(/i)/2 otherwise, 
whereas the critical field equals to he = JiA + J2{l + A)/2. 
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